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Abstract 

We provide detailed evidence for the claim that nonperturbative quantum grav- 
ity, defined through state sums of causal triangulated geometries, possesses a 
large-scale limit in which the dimension of spacetime is four and the dynamics of 
the volume of the universe behaves semiclassically. This is a first step in recon- 
structing the universe from a dynamical principle at the Planck scale, and at the 
same time provides a nontrivial consistency check of the method of causal dy- 
namical triangulations. A closer look at the quantum geometry reveals a number 
of highly nonclassical aspects, including a dynamical reduction of spacetime to 
two dimensions on short scales and a fractal structure of slices of constant time. 



1 



1 Introduction 



Nonperturbative quantum gravity can be defined as the quest for uncovering the 
true dynamical degrees of freedom of spacetime geometry at the very shortest 
scales. Because of the enormous quantum fluctuations predicted by the uncer- 
tainty relations, geometry near the Planck scale will be extremely rugged and 
nonclassical. Although different approaches to quantizing gravity do not agree 
on the precise nature of these fundamental excitations, or on how they can be 
determined, most of the popular formulations 1 agree that they are neither the 
smooth metrics 57^(2) (or equivalent classical field variables) of general relativity 
nor straightforward quantum analogues thereof. In such scenarios, one expects 
the metric to re-emerge as an appropriate description of spacetime geometry only 
at larger scales. 

Giving up the spacetime metric at the Planck scale does not mean discarding 
geometry altogether, since geometric properties such as the presence of a distance 
function pertain to much more general structures than differential manifolds with 
smooth metric assignments. One could hope that quantum gravity, when formu- 
lated in terms of the conjectured new microscopic degrees of freedom, was better 
behaved in the ultraviolet regime than standard perturbative approaches based 
on the spacetime metric. At the same time, it becomes a nontrivial test for 
such nonperturbative theories of quantum gravity whether they can reproduce 
the correct classical limit at sufficiently large scales. However, using this as a 
consistency check to discriminate between good and bad candidate theories is in 
practice complicated by the fact that some explicit results about the quantum 
dynamics of the proposed quantum gravity theory must be known. For example, 
this is not yet the case in loop quantum gravity jlj E] or in four-dimensional 
spin foam models for gravity [Zj. 

In the method of Causal Dynamical Triangulations one tries to construct a 
theory of quantum gravity as a suitable continuum limit of a superposition of 
spacetime geometries [HI HH ED] • In close analogy with Feynman's famous path 
integral for the nonrelativistic particle, one works with an intermediate regulariza- 
tion in which the geometries are piecewise flat 2 . The primary object of interest in 
this approach is the propagator between two boundary configurations (in the form 
of an initial and final spatial geometry), which contains the complete dynamical 
information about the quantum theory. Because of the calculational complexity 

exception are renormalization group approaches which look for nontrivial ultraviolet 
fixed points in a metric formulation P3E1E]. 

2 These are the analogues of the piecewise straight paths of Feynman's approach. How- 
ever, note that the geometric configurations of the quantum-gravitational path integral are not 
imbedded into a higher-dimensional space, and therefore their geometric properties such as 
piecewise flatness are intrinsic, unlike in the case of the particle paths. 
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of the full, nonperturbative sum over geometries (the "path integral"), an analyti- 
cal evaluation is at this stage out of reach. Nevertheless, powerful computational 
tools, developed in Euclidean quantum gravity [TTJ H21 HS1 UH EES 123 HZ] and 
other theories of random geometry (see ^H] for a review), can be brought to bear 
on the problem. 

This paper describes in detail how Monte Carlo simulations have been used 
to extract information about the quantum theory, and in particular, the geome- 
try of the quantum ground state 3 dynamically generated by superposing causal 
triangulations. It follows the announcement of several key results in this ap- 
proach to quantum gravity, first, a "quantum derivation" of the fact that space- 
time is macroscopically four-dimensional second, a demonstration that the 
large-scale dynamics of the spatial volume of the universe (the so-called "scale 
factor") observed in causal dynamical triangulations can be described by an ef- 
fective action closely related to standard quantum cosmology [20], and third, 
the discovery that in the limit of short distances, spacetime becomes effectively 
two-dimensional, indicating the presence of a dynamically generated ultraviolet 
cutoff [21] . They do not provide conclusive proof that our construction from first 
principles does lead to a viable theory of quantum gravity, but go some way in es- 
tablishing the existence of a physically meaningful classical limit. In addition, our 
detailed geometric analysis throws some light on the highly nonclassical micro- 
scopic properties of the spacetimes dominating the gravitational path integral, 
which after all are higher-dimensional analogues of the nowhere-differentiable 
paths that provide the support of Feynman's path integral. 

Our presentation will proceed as follows. After recollecting some basic prop- 
erties of the geometric set-up of the regularized path integral over causal trian- 
gulations in Sec. [21 we explain in Sec. [3] how the path integral is implemented 
numerically. The statistical model underlying our construction possesses three 
phases, which are distinguished by the geometric nature of their ground states, 
as described in Sec. [3] From this point on, we concentrate on the phase which 
yields a ground state which is extended in time and space. In Sec. 03 we provide 
the detailed numerical evidence for why this state is four-dimensional at large 
scales and why its shape can be described by a minisuperspace action. We then 
embark on a detailed investigation of substructures of the spacetime universe, 
namely, spatial slices (of fixed integer time r) in Sec. Eland in Sec. [7| "thick spa- 
tial slices", which are geometries of finite time extension At = 1. We measure 
both their Hausdorff and spectral dimensions, the latter from a discrete diffu- 
sion process. Interestingly, the geometry of the "thin" spatial slices shares some 
characteristics with a class of branched polymers which has already appeared in 

3 Here and in the following, by "ground state" we will always mean the state selected by 
Monte Carlo simulations, performed under the constraint that the volume of spacetime is (ap- 
proximately) kept fixed, a constraint we have to impose for simulation-technical reasons. 
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diverse as noncritical string theory, econophysics and network theory. A 
summary of our results and an outlook are presented in Sec. |H1 

2 Geometric set-up and simplicial action 

We briefly remind the reader of the structure of the spacetime geometries which 
contribute to the regularized path integral [10] . The fundamental building blocks 
are two types of four-simplices "cut out of four-dimensional flat Minkowski space. 
A flat four-simplex of type (4,1) has six space-like edges of length-squared /s pace = 
a 2 and four time-like edges of length-squared l 2 ime = —aa 2 . Similarly, a flat 
four-simplex of type (3,2) has four space-like edges of length-squared Zs pace = 
a 2 and six time-like ones of length-squared l 2 ime = —aa 2 . Each spacetime has 
a well-defined global foliation by a "proper time" r. 4 Fig. Q illustrates how 
the simplicial building blocks are situated between consecutive spatial slices of 
constant integer r. The (4,l)-simplex consists of a spatial tetrahedron at time 
r whose four vertices are connected by four time-like edges to the single vertex 
at time r + 1. Similarly, the (3,2)-simplex consists of a spatial triangle at time 
r whose vertices are connected to the two vertices at r + 1 which form the end 
points of the fourth spatial edge of the simplex. 

The way in which these building blocks, together with their counterparts with 
inversed time orientation, can form a spacetime slice of height At = 1 is by gluing 
them pairwise along a shared time-like boundary tetrahedron. In order to get a 
spacetime of time duration r = t, one glues together t such slices in succession, 
which is of course only possible whenever their spatial boundary geometries at 
integer-r match pairwise. We require the gluings to be such that the resulting 
discretized spacetime is a simplicial manifold. The spacetime topology is fixed to 
[0, 1] x where for our investigations we have chosen as constant-r surfaces 
the three-spheres ^T, = S 3 . For convenience in the computer simulations, we also 
sometimes identify time periodically, leading to a spacetime topology S 1 x S 3 . 

From the edge length assignments one can compute volumes and interior an- 
gles of the simplices which are needed because they appear in the Regge form 
of the Einstein action, from which the weight of each triangulated spacetime T 
in the sum over geometries will be computed (see ^H] for calculational details). 
The gravitational action becomes then a simple function of a number of positive 
integers, which count the overall numbers of simplices or subsimplices of a cer- 
tain type. For completeness, we give the explicit form of the action after analytic 

4 We choose the name "proper time" because the a priori integer-valued variable r can be 
extended continuously to the entire spacetime [55]. Within each simplex one can choose a 
coordinate system whose time variable coincides with r and also with the proper time of a 
freely falling observer in the simplex. Note that our setting here is that of piecewise flat and 
not of smooth geometries. 



4 



(4,1) 



(3,2) 



Figure 1: The two fundamental building blocks of causal dynamically triangulated 
gravity. The flat four-simplex of type (4,1) on the left has four of its vertices at time r 
and one at time t+1, and analogously for the (3,2)-simplex on the right. The "gap" 
between two consecutive spatial slices of constant integer time is filled by copies of 
these simplicial building blocks and their time-reversed counterparts, the (1,4)- and 
the (2,3)-simplices. 



continuation to the Euclidean sector, as a function of the number N of vertices 
(zero-simplices), and the numbers N^ 3 ' 2 ^ and AQ 4 ' 1 of four-simplices of types (3,2) 
and (4,1) of either time orientation, 
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In this expression, \ is the Euler characteristic of the piecewise flat four-manifold 
and d = —a denotes the positive ratio between the two types of squared edge 
lengths after the Euclideanization. In order to satisfy the triangle inequalities, we 
need d > 7/12 ^Uj- For simplicity, we have assumed that the manifold is compact 
without boundaries. In the presence of boundaries, appropriate boundary terms 
must be added to the action. 

For the simulations, a convenient alternative parametrization of the action is 
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given by 

Se = -(k + 6A)iV + k 4 (N^ + iVf l2) ) + A(2N[ 4 ' 1) + ivf' 2) ), (2) 

where the functional dependence of the and A on the bare inverse Newton 
constant k^ b \ the bare cosmological constant A® and a can be computed from 
(fUl. We have dropped the constant term proportional to Xi because it will be 
irrelevant for the quantum dynamics. Note that A = corresponds to 5 = 1, and 
A is therefore a measure of the asymmetry between the lengths of the spatial 
and timelike edges of the simplicial geometry As we shall see, the physically 
interesting region of the phase diagram of four- dimensional causal dynamical 
triangulations has positive A, which translates into values of a smaller than 1. 

3 Numerical implementation 

We have investigated the infinite-volume limit of the ensemble of causal triangu- 
lated four-dimensional geometries with the help of Monte Carlo simulations at 
finite four- volumes N 4 = ivj 4 ' 1 ^ + A^ 3 ' 2 ^ of up to 362.000 four-simplices. A simpli- 
cial geometry is stored in the computer as a set of lists, where the lists consist 
of dynamic sequences of labels for simplices of dimension n, < n < 4, together 
with their position and orientation with respect to the time direction. Additional 
list data include information about nearest neighbours, i.e. how the triangulation 
"hangs together" , and other discrete data (for example, how many four-simplices 
meet at a given edge) which help improve the acceptance rate of Monte Carlo 
moves. 

The simulation is set up to generate a random walk in the ensemble of causal 
geometries of a fixed time extension t. The local updating algorithm consists of 
a set of moves that change the geometry of the simplicial manifold locally, with- 
out altering its topological properties. These can be understood as a Lorentzian 
variant of (a simplified version of) the so-called Alexander moves [2H1 12H , i n 
the sense that they are compatible with the discrete time slicing of our causal 
geometries. For example, the subdivision of a four-simplex into five four-simplices 
by placing a new vertex at its centre is not allowed, because vertices can only be 
located at integer times r. Details of the local moves can be found in [10] . As 
usual, each suggested local change of triangulation is accepted or rejected accord- 
ing to certain probabilities depending on the change in the action and the local 
geometry. (Note that a move will always be rejected if the resulting triangula- 
tion violates the simplicial manifold property.) The moves are called in random 
order, with probabilities chosen in such a way as to ensure that the numbers of 
actually performed moves of each type are approximately equal. We attained a 
rather high average acceptance rate of about 12.5%, which was made possible by 
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keeping track of sub-simplex structures relevant for the performance of certain 
moves, as mentioned earlier. The attempted moves were performed in units of 
10 6 (one "sweep"), with the time for one sweep approximately independent of 
the system volume in the region k,q « 2.2 and A « 0.4 . . . 0.6 where most of our 
measurements were taken. 

In order to obtain an efficient sampling of geometries of large four- volume (the 
larger the four- volume, the better the computer-generated geometries describe the 
continuum limit), we have performed the simulations at (approximately) constant 
four-volume. This means that instead of the standard gravitational path integral 
(after Euclideanization) 

Z E (X,G) = J V[g] e- s ^\ S E [g] = -1 J d 4 x v^det^ (R - 2A) (3) 

we simulate a discretized version of the path integral at fixed four- volume V4, 

Z E (V 4 ,G) = J V[g}e-~ S ^H(J d A x^~g -V,), S E [g] =~ J dPxy/detg R, 

(4) 

which is related to (|3*|) by a Laplace transformation, 

POO 

Z E {\, G) = / dV 4 e-^Z E (V,, G). (5) 
Jo 

However, the local moves do not in general preserve the numbers N^' 1 ' and N^' 2 ' 
of four-simplices, or indeed their sum. We deal with this in a standard way which 
was developed for dynamically triangulated models in dimensions three and four 
E] to ensure that the system volume is peaked at a prescribed value, with a 
well-defined range of fluctuations. Adapting it to the current situation of causal 
triangulations with nonvanishing asymmetry A, we implement an approximate 
four- volume constraint by adding a term 

<J5 , = e|JVi 4 ' 1) -JV 4 |, (6) 

to the Euclidean action, with typical values of e lying in the range of 0.01 to 0.02, 
except during thermalization where we set e = 0.05. The reason for fixing A"] 4 ' 1 "* 
instead of A" 4 = A"| 4>1) + A"| 3,2) in eq. (0) is mere technical convenience. We have 
checked in the phase space region relevant to four- dimensional quantum gravity 
(phase C, see below) that for A^ 4 ' 1 ^ fixed according to (JUJ), the number A^ 3 '"' 
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of four-simplices of type (3,2) is likewise very sharply peaked, see Fig. The 
"four-volumes" A^ 4 and the corresponding "true" discrete four-volumes A" 4 used 
in the simulations are listed in Table ^ In order to stabilize the total volume 
after thermalization, k 4 has to be fine-tuned to its pseudo-critical value (which 
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"four-volume" N 4 = ivj 4,1) 


10 


20 


40 


80 


160 


actual four- volume N4 = +N^' 2 ^ 


22.25 


45.5 


91 


181 


362 



Table 1: Translation table between the two types of discrete four- volume, N4 and 
N±, at which the numerical simulations reported in this article were performed, 
in units of 1000 building blocks (phase C only). 

depends weakly on the volume) with accuracy smaller than e, in practice to about 
0.2e. The measurements reported in this paper were taken at N± = 10, 20, 40, 80 
and 160k, and the runs were performed on individual PCs or a PC farm for the 
smaller systems and a cluster of work stations for the larger systems. 

Before measurements can be performed, one needs a well thermalized con- 
figuration of a given volume. In order to double-check the quality of the ther- 
malization, we used two different methods to produce starting configurations for 
the measurement runs. In the first method, we evolved from an initial mini- 
mal four-dimensional triangulation of prescribed topology and of a given time 
extension t, obtained by repeated gluing of a particular triangulated space-time 
slice of At — 1 and topology [0, 1] x S* 3 , which consists of 30 four-simplices. The 
spatial in- and out-geometries of the slice are minimal spheres S 3 , made of five 
tetrahedra. The two types of spatial boundary conditions used are (i) periodic 
identification of the geometries at initial and final integer times, and (ii) free 
boundary conditions, where all vertices contained in the initial slice at time To 
are connected by timelike edges to a single vertex at time To— 1, and similarly for 
the vertices contained in the final spatial slice. From this initial configuration, 
the geometry evolves to its target volume N4 as specified through 5S. During 
the evolution the volume-volume correlator (see eq. (JJJ) below) changes from a 
very broad distribution to the stable narrow shape seen in Fig. [7| The number 
of sweeps to reach the thermalized configuration changes linearly with N4 and 
ranged from 10 5 to 10 8 sweeps for the largest volumes, the latter of which took 
several weeks on a work station. 

In the second method, we started instead from a thermalized configuration of 
smaller volume, which we let evolve towards the target volume. In this case the 
final volume- volume distribution is reached from a narrower distribution, namely, 
that of the smaller volume. During thermalization, this width grows very slowly. 
The timing of the entire process is similar to that of the first method. 
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Figure 2: Unnormalized distribution of the number AQ ' of four-dimensional simplices 



of type (3,2), at fixed numbers of four-simplices of type (4,1), for iV 
160k (left to right) at k = 2.2 and A = 0.6. 



(4,1) 



40, 80 and 



4 Phase structure of the model 

Like always in dynamically triangulated models, the bare cosmological constant 
(equivalently, k 4 ) is tuned to its (pseudo-)critical value in the simulations, 
tantamount to approaching the infinite- volume limit. Depending on the values 
of the two remaining parameters k,q and A in the discretized action ©, we have 
identified three different phases, A, B and C, mutually separated by lines of first- 
order transitions. We have not yet measured their location in detail; Fig. |3] is a 
qualitative sketch of the phase diagram. We describe each of the three phases in 
turn: 

(A) This phase prevails for sufficiently large k (recall k q is proportional to the 
bare inverse Newton's constant k^ b '). When plotting the volume of the 
spatial slices r = const as a function of r, we observe an irregular sequence 
of maxima and minima, where the minimal size is of the order of the cutoff, 
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A 



B 



A 



Figure 3: A qualitative sketch of the phase diagram of four-dimensional causal dy- 
namical triangulations as function of the gravitational coupling kq and the asymmetry 
parameter A. 



and the sizes of the maxima vary, see Fig. The time intervals during 
which the spacetime has a macroscopic spatial extension are small and of 
the order of At = 3. 

(B) This phase occurs for sufficiently small kq and for small asymmetry A, 
including A = 0. In it, spacetime undergoes a "spontaneous dimensional 
reduction" in the sense that all four-simplices are concentrated in a slice of 
minimal time extension Ar = 2, and the three- volume N 3 (t) remains close 
to its kinematic minimum everywhere else (Fig. [SJ. 

We believe that these two phases can be understood essentially from the dynamics 
of the (Euclidean) geometries of the slices of constant r alone. The spatial slices 
in phase A are presumably realizations of branched polymers which dominate 
the path integral of three-dimensional Euclidean quantum gravity for large kq, 
whereas the extended part of the geometry in phase B corresponds to a Euclidean 
geometry in the so-called crumpled phase, which is characterized by the presence 
of very few vertices with extremely high coordination number j2Hl El [13 Ell • 
We will not elaborate here on our detailed measurements which quantify these 
statements, because we view phases A and B as lattice artifacts. However, there is 
yet another phase where the dynamics is not reduced to just the spatial directions 
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Figure 4: Monte Carlo snapshot of a typical universe in phase A (kq = 5.0, A = 0), 
of discrete volume N& = 45.5k and total time extent (vertical direction) t = 20. In 
this and the following two figures, the circumference at integer proper time r is chosen 
proportional to the spatial three- volume V%{t). The surface represents an interpolation 
between adjacent spatial volumes, without capturing the actual 4d connectivity between 
neighbouring spatial slices. 

but is genuinely four- dimensional. 

(C) This phase occurs for sufficiently small Kq and nonvanishing asymmetry A. 
In this phase, where A > and a < 1 (5 was denned below eq. ((TJ)), there 
seems to be a sufficiently strong coupling between successive spatial slices 
to induce a change in the spatial structure itself. We will discuss this new 
geometrical structure in detail in Sec. El 

Of course, the presence of the asymmetry parameter A is only possible be- 
cause of the Lorentzian character of the model, and the resulting distinction 
between spatial and time directions. Since we implemented this distinction 
in the discretization by setting up a particular proper-time foliation and 
choosing arbitrary discrete length and time scales, it should not be surpris- 
ing that in order to make contact with continuum physics, the ratio of these 
two scales cannot be chosen completely arbitrarily. This merely implies a 
restriction on the way the discretization is set up, and does not lead to any 
new parameters in the continuum theory. We have checked that the physi- 
cal results derived within phase C do not depend on the numerical value of 
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Figure 5: Monte Carlo snapshot of a typical universe in phase B (kq = 1.6), of discrete 
volume ^4 = 22.25k and total time extent t = 20. The entire universe has collapsed into 
a slice of minimal time extension. 

A, within the allowed range. As we will see later, any such choice leads to 
the large-scale geometry of a homogeneous and isotropic universe. 

The hallmark of this phase is the presence of a stable extended four-geometry, 
as first reported in jTH] , and illustrated in Fig. |U1 Not only does the macro- 
scopic dimensionality of spacetime emerge with the correct classical value, 
but, equally remarkably, the global shape of the universe has been shown 
to be related to a simple minisuperspace action similar to those used in 
standard quantum cosmology [2D]- Details of the computer measurements 
in this phase are the subject of Sees. 030 below. 

Looking at the phase diagram of Fig. El it may at first be surprising that one 
does not obtain a continuum limit when the bare gravitational coupling constant 
(or Newton's constant) vanishes, G^ b ' ~ — ► 0. However, one should bear in 
mind that the method of dynamical triangulation is a way of defining the path 
integral, rather than a regularization method which has nice smooth continuum 
manifolds as its dominant configurations as — > 0. 

A helpful analogy is the example of the (Euclidean) path integral for the 
relativistic particle, whose propagator G(x, y) can be represented as a sum over all 
random paths from x to y, whose action S(P) = m^ b 'L(P) is given by the length L 
of the path P in target space, multiplied by a dimensionless bare mass parameter 
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m^ h \ If this path integral is regularized on a hypercubic lattice in R d , the possible 
paths are by definition along the lattice links. Now consider R 2 and a path from 
(0, 0) to the lattice point (x, y) = (a ■ n,a ■ n). No matter how small the lattice 
spacing a and how large the number n of lattice steps, the lattice paths between 
the two points will never approach smooth continuum paths and their lengths 
will always be larger than \[2\J x 2 + y 2 . To define a continuum limit for the 
lattice random walk one looks for critical points of the corresponding statistical 
system and studies the long-distance properties of suitable observables, in this 
case, of the propagator. The critical point is determined by the value of mS b \ 
but this bare coupling constant is usually not related to the physical parameter 
(given here by the mass term in the continuum propagator) in a straightforward 
way (see JHIEZI) for details). In the case of dynamical triangulations we follow 
a similar strategy. We do not a priori associate a continuum limit with certain 
values of the bare coupling constants. Rather, we look for points or regions in 
coupling-constant space where interesting scaling limits may exist. In the case 
of four-dimensional gravity, it is natural to look for regions where some desirable 
features of macroscopic spacetime are reproduced, and this is precisely the case 
for the phase-space region (C). 

Let us also comment on how the continuum limit is taken in our approach. 
According to the standard, Wilsonian understanding of the relation between the 
theory of critical phenomena in statistical mechanics and quantum field theory, in 
order to obtain a continuum quantum field theory from an underlying regularized 
statistical system one needs to fine-tune a few coupling constants (correspond- 
ing to the relevant directions in coupling constant space) to reach the infinite- 
dimensional critical surface where a suitably defined correlation length is infinite. 
Directions along the critical surface correspond to the infinitely many irrelevant 
(nonrenormalizable) directions in coupling constant space. Weinberg conjectured 
that this picture, which is usually believed to define ordinary renormalizable 
quantum field theories, may generalize to quantum gravity, in the sense that it 
may have an associated, nonperturbatively defined critical surface |28j . Several 
groups have since tried to substantiate this idea using either a (2 + e)-expansion 
[29J or refined renormalization group techniques jSj- Various lattice approaches 
to quantum gravity like quantum Regge calculus |3U] , Euclidean dynamical trian- 
gulations [TTJ EH C2] and the causal dynamical triangulation of quantum gravity 
[TU1 1201 UH] used by us presently also endorse this picture. 

In the latter, Lorentzian approach, the bare cosmological constant must be 
fine-tuned to obtain an infinite spacetime volume. Formally, we are not perform- 
ing an analogous fine-tuning of the bare gravitational coupling constant in phase 
C, which is nonstandard from a Wilsonian point of view. There are indeed lat- 
tice theories where no fine-tuning is needed, and where an infinite- volume limit 
is sufficient to identify a continuum limit, for example, the Coulomb phase of 
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[/(l)-lattice gauge theory. However, the situation in gravity could be more in- 
volved, and subtleties may be buried in the identification of more complicated 
observables than the global ones we have considered so far. For a related dis- 
cussion in three-dimensional causal dynamically triangulated quantum gravity, 
we refer to pffi 133] . What is important from our present point of view is the 
fact that the universes emerging from phase C satisfy the first nontrivial tests of 
being viable candidates for quantum universes, as we will now turn to describe. 
Encouraged by these results, we are in the process of developing new and more 
refined methods to probe the geometry of the model further, and which eventu- 
ally should allow us to test aspects related to its local "transverse" degrees of 
freedom, the gravitons. We invite and challenge our readers to find such tests in 
a truly background-independent formalism of quantum gravity. 

5 Measuring geometry 

Having generated a ground state of quantum geometry nonperturbatively, we 
would like to understand its geometric properties, in the sense of expectation 
values, which can be done in various ways. We will currently concentrate on the 
purely geometric observables, leaving the coupling to test particles and matter 
fields to a later investigation. We will proceed by first determining a number 
of "rough" properties of the quantum geometry such as its dimensionality and 
its global scaling properties. These properties of quantum spacetime play an 
important role, because for any viable candidate theory of quantum gravity, they 
must be shown to reproduce the correct classical limit at sufficiently large scales, 
namely, geometry as described by Einstein's general theory of relativity. 

As we have emphasized in ^J] , i n nonperturbative models of quantum gravity, 
where by definition geometry fluctuates very strongly at short scales, there is 
absolutely no guarantee that any sensible classical limit is obtained even for very 
basic quantities like the dimension of spacetime. Indeed, previous formulations 
which can probe the dimensionality of geometry numerically through the scaling 
behaviour of suitable observables, for example, dynamically triangulated models 
of Euclidean quantum gravity, have not come up with evidence for the value 4. 
Instead such models seem to run generically into highly degenerate geometries 
of one type or other ^21 EU El EE] > which at large scales have nothing to do 
with four-dimensional spacetime as we know it. From this it seems likely that the 
condition of obtaining the correct classical limit imposes strong constraints on 
any attempt to quantize gravity nonperturbatively and background-independently. 
Of course, proposals of quantizing gravity are often not sufficiently developed 
to decide whether they can reproduce classical geometry from first principles, 
without putting in a preferred background structure along the way. 
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In the following, we will present a detailed discussion of a number of geometric 
properties of the extended ground state of phase C above. Several main results 
were already announced in our previous publications [T^l 1201 EE] • The picture 
emerging so far is that of a geometry whose large-scale properties are classical, 
but whose detailed microscopic structure is highly complicated and nonclassi- 
cal. We will proceed by first presenting the evidence for the macroscopic four- 
dimensionality of spacetime, and then describing the insights into the geometric 
structure obtained from measuring quantities associated with slices of constant 
time. 




Figure 6: Monte Carlo snapshot of a typical universe in phase C («o=2.2) of discrete 
volume ^4 = 91.1k and total time extent t = 40. 



5.1 The evidence for spacetime dimension four: 
Global scaling 

As illustrated by Fig. El the total four-volume is distributed over the t time 
steps (we are working with a periodic identification in the time direction) in a 
characteristic way. At a given discrete four-volume, for sufficiently large t the 
configuration separates out into a thin "stalk" whose four-volume at each time 
step fluctuates little and is close to the kinematically allowed minimum for the 
topology [0, 1] x S 3 , and an extended, blob-shaped structure, the "universe". The 
average volume s per time step At = 1 of the stalk is practically independent of the 
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total volume of the system, with s ~ 25.5 for the range of coupling constants we 
used. As explained in Sec. 01 we measure the spacetime volume by iV 4 , that is, in 
units of four-simplices of type (4,1). Since every spatial tetrahedron is contained 
in exactly two (4,l)-simplices, the average three- volume A^(r) at integer time r 
in the stalk is therefore given by s/2. 

A suitable observable for studying the large-scale structure of the universe is 
the volume-volume correlator 5 

c«. (S) , Mows)) = ± W)-./|)Wr + f )-./2 )> | (7) 

which satisfies the normalization condition 

t-i 

<5=0 

where a periodic identification of the time variable has been assumed. Because 
of the subtraction of the cutoff volume s/2, the volume- volume correlator is only 
sensitive to the extended part of the geometry, and measures the distribution of 
the "effective" four-volume 

Nf = A> 4 - ts (9) 

contained in it. We have measured the distribution (JJJ) for system sizes N^ = 10, 
20, 40, 80 and 160k, at k = 2.2, A = 0.6 and t = 80. We expect the distribution 
to scale as a function of N4, with an underlying universal distribution 

(Nfy/^c^uNfy/^x) (io) 



C AT , [X 



JV 4 W ■- V 1V 4 ) N4, 

depending on the rescaled time variable 6 

x = _ !? . (11) 

The value of the scaling dimension Dh is then determined by the method of 
"best overlap". In it, we compare each measured distribution with that of the 
40k system. In order to provide a comparison point for each value of x and 
to quantify the overlap, we first performed a spline interpolation of the discrete 
distribution at 40k. Fig. [7| shows a combined plot of the rescaled distributions, 
using Dh = 4 as a scaling dimension, and with a plot range symmetrized around 
x — 0. The overlap is nearly perfect. As can be read off Fig. [HI which shows the 
average error of the overlap as a function of Dh, the best overlap corresponds to 
Dh ~ 3.8 ± 0.25. We take this outcome as strong evidence that our spacetimes 

5 Note that in the analogous relation (6) in 19 we mistakenly omitted the ensemble averages 
on the right-hand side. 

6 We are employing standard finite-size scaling methods familiar from the theory of critical 
phenomena (see [HE] for a textbook exposition), adjusted to our problem. 
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Figure 7: The scaling of the volume-volume distribution, as function of the rescaled 
time variable x = S/(N^) 1 ^. Data points come from system sizes N^ = 10, 20, 40, 80 
and 160k at k = 2.2, A = 0.6 and £ = 80. 

are genuinely four- dimensional at large scales. 7 

5.2 The evidence for spacetime dimension four: 
Spectral dimension 

Another way of obtaining an effective dimension of the nonperturbative ground 
state, its so-called spectral dimension D$, is from studying a diffusion process on 
the underlying geometric ensemble. On a d- dimensional manifold with a fixed, 
smooth Riemannian metric g a b{£), the diffusion equation has the form 

^K g (Uo^) = A g K g {Uo'^), (12) 

where a is a fictitious diffusion time, A g the Laplace operator of the metric g a b(0 
and K g (£, £o; o) the probability density of diffusion from £ to £o in diffusion time 

7 Of course, we will never be able to show by computer simulations alone that the effective 
dimension is exactly equal to four. What we are assuming here is that this dimension is 
indeed an integer, because there are no classical theories describing the dynamics of geometries 
of dimensionality 3.9, say. Also, it is known from the measurement of Hausdorff dimensions 
in two-dimensional Euclidean quantum gravity (where it can be proven that the Hausdorff 
dimension is four (!) |3fil I37| ) that the computer simulations tend to slightly underestimate the 
actual values CHH ■ 



17 



0.002 



.0015 



.001 



.0005 



3 3.25 3.5 3.75 4 4.25 4.5 4.75 

D H 

Figure 8: Estimating the best scaling dimension Dh from optimal overlaps of rescaled 
volume- volume distributions. The plot variable is the average error of the overlap as 
function of Dh- 



a. We will consider diffusion processes which initially are peaked at some point 



2U£,£o;<7=0) 



1 



5 d (£-£o). 



For the special case of a flat Euclidean metric, we have 



e -<?(f.6>)/*T 



g a b(0=dab- 



For general curved spaces K g has the well-known asymptotic expansion 



e -4(^€o)/4c7 ~ 

Kgfa £o! o-) ~ 2^ a r (f , Co) ^ r 

r=0 



a 



(13) 



(14) 



(15) 



for small a, where d g (£, £o) denotes the geodesic distance between £ and £ - Note 
the appearance of the power a~ d ^ 2 in this relation, reflecting the dimension d of 
the manifold, just like in the formula ()14j) for flat space. This happens because 
small cr's correspond to short distances and for any given smooth metric short 
distances imply approximate flatness. 

A quantity that is easier to measure in numerical simulations is the average 
return probability P g {cr), which possesses an analogous expansion for small a, 



If 1 °° 



(16) 



where V is the spacetime volume V = J d d ^ w det g (f ) and the expansion coeffi- 
cients A r are given by 

A r = 1^ d^ydet^O a r (e,0- (17) 

For an infinite flat space, we have P g (a) = l/(47ra) d / 2 and thus can extract the 
dimension d by taking the logarithmic derivative, 

_ 2 dlogP» 
a logo- 
independent of a. For nonflat spaces and/or finite volume V, one can still use eq. 
(|18)1 to extract the dimension, but there will be corrections for sufficiently large 
o. For finite volume in particular, P g (cr) goes to 1/V for a ^> V 2 l d since the 
zero mode of the Laplacian — A g will dominate the diffusion in this region. For 
a given diffusion time a the behaviour of P g (cr) is determined by eigenvalues A n 
of — A g with A n < I/a - , and the contribution from higher eigenvalues is exponen- 
tially suppressed. Like in the flat case, where diffusion over a time a probes the 
geometry at a linear scale a/cx, large a corresponds to large distances away from 
the origin £ of the diffusion process, and small a to short distances. 

The construction above can be illustrated by the simplest example of diffusion 
in one dimension. The solution to the diffusion equation on the real axis is 

e"« 2 / 4CT 

Jjf(e,a) = -=, K(k,a) = e- k °, (19) 

where K(k,a) denotes the Fourier transform of K(£,a). The eigenvalues of the 
Laplace operator are of course just given by k 2 . In order to illustrate the finite- 
volume effect, let us compactify the real line to a circle of length L. The return 
probability is now given by 



where in the second step we have performed a Poisson resummation to highlight 
the cr -1 / 2 -dependence for small a. In line with the discussion above, the associated 
spectral dimension D s (u) is constant (and equal to one) up to a- values of the 
order L 2 /4ir 2 , and then goes to zero monotonically. 

In applying this set-up to four- dimensional quantum gravity in a path integral 
formulation, we are interested in measuring the expectation value of the return 
probability P g (<y). Since P g {o) defined in (fi7)J) is invariant under reparametriza- 
tions, it makes sense to take its quantum average over all geometries of a given 
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spacetime volume V4, 



1 



e- §B ^6{ / d 4 x V / dit^ - V 4 ) P a {a) 



(21) 



where the partition function Ze and the corresponding gravitational action Se 
were defined in Since the small-cr behaviour of -P g (cr) is the same for each 
smooth geometry, it might seem obvious that the same is true for their integral 
Py 4 (<r), but this need not be so. Firstly, the scale a in (f2~Tj) is held fixed, in- 
dependent of the geometry g ab , while the expansion (|16p contains reference to 
higher powers of the curvature of g a b- Secondly, one should keep in mind that a 
typical geometry which contributes to the path integral - although continuous - 
is unlikely to be smooth. This does not invalidate our treatment, since diffusion 
processes can be meaningfully defined on much more general objects than smooth 
manifolds. For example, the return probability for diffusion on fractal structures 
is well studied in statistical physics and takes the form 



where N is the "volume" associated with the fractal structure and Ds the so- 
called spectral dimension, which is not necessarily an integer. An example of 
fractal structures are branched polymers (which will play a role in Sec. |H1 below 
when we analyze the geometry of spatial slices), which generically have D$ = 4/3 
jJOlE]- Extensive numerical simulations OllSj have shown that in 2d quantum 
gravity the only effect of integrating over geometries is to replace the asymptotic 
expansion pfij) . which contains reference to powers of the curvature related to a 
specific metric, by the simpler form ([22)1 . 

Our next task is to define diffusion on the class of metric spaces under con- 
sideration, the piecewise linear structures defined by the causal triangulations T. 
We start from an initial probability distribution 



which vanishes everywhere except at a randomly chosen (4,1) -simplex iq, and 
define the diffusion process by the evolution rule 



These equations are the simplicial analogues of (fT3J) and (|T2*jl . with the triangu- 
lation (together with its Euclideanized edge-length assignments) playing the role 




(22) 




(23) 




(24) 
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of g a b, and k — > j denoting the five nearest neighbours of the four-simplex j. In 
this process, the total probability 

Y^K T (j,i ',*) = l (25) 

3 

is conserved. The return probability to a simplex i is then defined as Pr(io] cr) = 
Kt^qi io', o") and the quantum average as 

= Tl^T,^'^ IF E (26) 



y ' T N4 «oGTjv 4 



where T/v 4 denotes a triangulation with N4 four-simplices, and Se(Tn 4 ) and 
Ze(N±) are the obvious simplicial analogues of the continuum quantities (HJ) at 
fixed four-volume. Assuming that the return probability behaves according to 
(I22|) . with N = N^, we can extract the value of the fractal dimension D$ by mea- 
suring the logarithmic derivative as in (JTSj) above, as long as the diffusion time is 
not much larger than N^ Ds , 

D s (a) = —2 — ^ — — - + finite-size corrections. (27) 
a log a 

From the experience with numerical simulations of 2d Euclidean quantum 
gravity in terms of dynamical triangulations El HB] , we expect some irreg- 
ularities in the behaviour of the return probability for the smallest a, i.e. close 
to the cut-off scale. Typically, the behaviour of Pjv(c) for odd and even diffusion 
steps cr will be quite different for small a and merge only for a ~ 20 — 30. After 
the merger, the curve enters a long and stable regime where the right-hand side 
of (J2*7j) is independent of a, before finite-size effects start to dominate which force 
it to go to zero. (Also the 3d Euclidean hypersurfaces of constant time in our 
current set-up show a very similar behaviour, as can be seen in Fig. [TBI) 

The origin of the odd-even asymmetry can again be illustrated by the simple 
case of diffusion on a one-dimensional circle, whose solution is the discretized 
version of the solution (J20j) . In this case the asymmetry of the return probability 
between odd and even time steps is extreme: if we use the two-dimensional version 
of the simple evolution equation (|24|). we obtain 

for cr odd, 

Pl(°) =\ 1 f a \ (28) 
' for cr even, 



2 CT \(j/2 

as long as cr < L/2, where L is the discrete volume of the circle (i.e. the number 
of its edges). It is of course possible to eliminate this asymmetry by using an 
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Figure 9: The spectral dimension Dg{a) of causal dynamical triangulations as a func- 
tion of the diffusion time a, which is a direct measure of the distance scale probed. 
The measurements were taken at volumes A^ = 20k (bottom curve), 40k and 80k (top 
curve), and for kq = 2.2, A = 0.6 and £ = 80. 

"improved" discretized diffusion equation, but in the case of higher-dimensional 
random geometries like the ones used in 4d causal dynamical triangulations this 
is not really necessary. The random connectivity will eliminate the asymmetry 
after a number of steps as exemplified by Fig. 

The results of measuring the spacetime spectral dimension D$ were first re- 
ported in [21] . We work with system sizes of up to N± = 80k with kq = 2.2, A = 0.6 
and £ = 80. Since we are interested in the bulk properties of quantum spacetime 
and since the volume is not distributed evenly in the time direction (cf. Fig. HJ) , 
we always start the diffusion process from a four-simplex adjacent to the slice 
of maximal three-volume. When this is done the variations in the curve Ds(<j) 
for different generic choices of Tjy and starting simplices Iq are small, as long 
as io is located not too close to the stalk. The data curves presented in Fig. El 
were obtained by averaging over 400 different diffusion processes performed on 
independent configurations. Because following the diffusion out to a time a = 400 
is rather time-consuming and the variation in Pt^o', cr) as a function of T and 
io (with the initial conditions specified above) is small, we have not performed 
a more extensive average over configurations. We have omitted error bars from 
Fig. El to illustrate how the curves converge to a shape that represents Ds{cr) in 
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the infinite- volume limit, and which is given by the envelope of the data curves 
for finite volume. For the two lower volumes, N4 = 20k and N4 = 40k, there still 
are clear finite-volume effects for large diffusion times a. 

By contrast, the top curve - corresponding to the maximal volume A/4 = 80k 
- continues to rise for increasing a, which makes it plausible that we can ignore 
any finite-size effects and that it is a good representative of the infinite-volume 
limit in the cr-range considered. 8 We will therefore concentrate on analysing the 
shape of this curve, which is presented separately in Fig. now with error bars 
included. (More precisely, the two outer curves represent the envelopes to the 
tops and bottoms of the error bars.) The error grows linearly with a, due to the 
occurrence of the logcr in (|27j) . 

The remarkable feature of the curve Ds(cr) is its slow approach to the asymp- 
totic value of Ds(cr) for large a. As emphasized in [21], this type of behaviour has 
never been observed previously in systems of random geometry, and again under- 
lines that causal dynamical triangulations in four dimensions behave qualitatively 
differently, and that the quantum geometry produced is in general richer. The 
new phenomenon we observe here is a scale dependence of the spectral dimension, 
which has emerged dynamically. This is to be contrasted with fractal structures 
which show a self-similar behaviour at all scales, and which lead to -Ds(cx)-curves 
like that shown in Fig. 

As explained in [21], the best three-parameter fit which asymptotically ap- 
proaches a constant is of the form 

b 119 
D s (a) = a-—- = 4.02 - ——. (29) 
a + c 54 + a 

The constants a, b and c have been determined by using the full data range 
a G [40,400] and the curve shape agrees well with the measurements, as can be 
seen from Fig. Integrating (|29|) we obtain 

m ~ ^ + 1 c/g)V2e , (30) 

from which we deduce the limiting cases 

o~ a l 2 for large a, 
FV)H ' (31) 

a -(a-b/c)/2 for gmall a _ 



8 In both Figs. El and we have only plotted the region where the curves for odd and even 
a coincide, in order to exclude short-distance lattice artifacts. The two curves merge at about 
a = 40. Since the diffusion distance grows as y/a, a return diffusion time of 40 corresponds to 
just a few steps away from the initial four-simplex. 
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Figure 10: The spectral dimension D$ of the universe as function of the diffusion time 
a, measured for ko = 2.2, A = 0.6 and £ = 80, and a spacetime volume ^4 = 181k. The 
averaged measurements lie along the central curve, together with a superimposed best 
fit Ds(cr) = 4.02 — 119/(54+0"). The two outer curves represent error bars. 

We conclude that the quantum geometry generated by causal dynamical trian- 
gulations has a scale-dependent spectral dimension which increases continuously 
from a — b/c to a with increasing distance. Substituting the values for a, b and c 
obtained from the fit (|29jh and taking into account their variation as we vary the 
<7-range [cr min , cr max ] and use different weightings for the errors, we obtained the 
asymptotic values [2T] 



for the "short-distance spectral dimension". A dynamically generated scale- 
dependent dimension with this behaviour is truly exciting news, because it sig- 
nals the existence of an effective ultraviolet cut-off for theories of gravity with 
and without matter coupling, brought about by the (highly non-perturbative) be- 
haviour of the quantum-geometric degrees of freedom on the very smallest scale. 

5.3 The effective action 

In order to get a more detailed understanding of the geometric properties of our 
universe, we next attempt to reproduce the overall shape of the universe of Fig. 




(32) 



for the "long-distance spectral dimension" and 




(33) 
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|B] (in other words, the behaviour of its three-volume V3 as a function of time) 
from an effective action for the single variable V^t), valid at sufficiently large 
scales V3 ^> 1. This will enable us to compare our results with those from so- 
called minisuperspace models. The simplest such models also take the form of 
dynamical systems of a single variable, namely, the scale factor a(r), describing 
the linear spatial extension of the universe and related to the three-volume by 

V 3 (t) oc a 3 (r). (34) 

Despite the resemblance, the status of our derivation is very different from that 
of a minisuperspace model. In the latter, one works from the outset in a dras- 
tically truncated framework, where the classical dynamics of the entire universe 
is described by the variable a(r) instead of the full metric data gij(r, x). For 
closed universes, and in (Euclidean) proper-time coordinates, this corresponds to 
considering line elements of the form 

d s 2 = dr 2 + a 2 (r)dnl, (35) 

where dfl 2 denotes the metric on the three-sphere. The corresponding Einstein- 
Hilbert action 9 is 

5— = I J dr (-a(r) ' - a{r) + Aa 3 (r) j , (36) 

including a cosmological constant A. These minisuperspace models are often taken 
as starting point for a path integral quantization, leading to quantum cosmologies 
like those used by Hartle and Hawking in their semiclassical evaluation of the 
wave function of the universe US] (see also 0HI 113 HH1 113 E01 EU E21 EHJ for 
related work). The hope is that such a truncated version of quantum gravity still 
describes some features of the full theory correctly. 

By contrast, the central feature of dynamical triangulations is a bona fide path 
integration over all geometric degrees of freedom (which in particular will have 
a nontrivial spatial dependence), without appealing to any symmetry reduction, 
either classical or quantum-mechanical. In the final result, one may then decide to 
keep track only of the three- volume (equivalently, the scale factor) as a dynamical 
variable and try to extract its effective large-scale behaviour. It is an interesting 
question whether the two approaches (that is, truncating the degrees of freedom 
before or after the quantization) are equivalent. In view of the severe difficulties 
in making quantum-cosmological path integrals well-defined and convergent [ST] , 
one might hope that they are not. This is indeed the conclusion from the findings 

9 We work with Euclidean signature in order to facilitate comparison with the simulation 
data, which also describe Euclidean, Wick-rotated geometries. 
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in causal dynamical triangulations we are about to describe. As we will see, the 
collective effect of the full quantum-gravitational degrees of freedom changes the 
dynamical description in terms of the single scale factor in a drastic way. 
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Figure 11: Extracting the scaling dimension D 2 by finite-size scaling, in order to 
determine the kinetic term in the effective action of the scale factor. 

Our first step is to determine the kinetic term in the effective action of the 
dynamically triangulated model, and compare it to the term —aa 2 of the mini- 
superspace model ()36j) . We obtain this term by measuring the distribution of 
differences of adjacent spatial volumes 

z JN 3{ r + t>-VN 3{T)[ w(T)+iV3(T+1) . (37) 

for different values of the discrete three- volume N3. Again, the information we 
extract is two-fold: we first determine the exponent D 2 from best overlap, and 
subsequently the overall shape of the universal distribution. In order to maximize 
the sampled three-volumes, we employ the free boundary conditions described in 
Sec. EJabove. That is, all configurations extend over three time steps, At = 3, with 
only a single vertex each at the initial time r = and the final time t = 3, and the 
relevant volume difference given by (|37p. with r = l. From taking measurements 
at system volumes N4 = 10, 20, 40, 80 and 160k, we extracted a best value for 
D 2 by minimizing the failure of the curves to overlap for different volumes, using 
the same method described earlier for the determination of the scaling dimension 
Dh ■ The result is illustrated by the curve of Fig. which has its minimum at 
about 2.12, which we interpret as evidence for D 2 = 2. The overlapping curves 
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for the volume distributions Pj^ 3 (z) at D 2 = 2 are depicted in Fig. H21 Apart from 
^-values near zero, the overlap is clearly good. Moreover, the curve is fitted very 
well by a Gaussian e~ cz , with a positive constant c independent of N 3 . 
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Figure 12: The distribution Pn 3 (z) of rescaled volume differences z between adjacent 
spatial slices, for three-volumes -/V3 = 20, 40, 80 and 160k. 

From estimating the entropy of spatial geometries, that is, the number J\f(N 3 ) 
of three-dimensional configurations at given three- volume N 3 , one expects jUj 
corrections of the form — c'iVJ, with < a < 1, d > 0, to the exponent —cz 2 = 
— c (AA^ 3 /Ar) 2 /A^3 in the distribution P^ 3 (z). This leads us to the following form 
for the Euclidean effective action, valid for large three-volume N 3 (t), 

*S « t (iirj (^) 2 + ^"M - A ^>) • < 38 > 

with < a < 1, and ci, C2 > 0. In (jSHJ), A is a Lagrange multiplier to be 
determined such that 

t 

J2n 3 (t)=N 4 . (39) 

r=0 

Unfortunately it is impossible to measure the precise form of the correction term 
C2N3 directly in a reliable way. However, it is clear from a general scaling of 
the above action that the only possibility for obtaining the observed scaling law, 
expressed in terms of the variable t/N\ , is by setting a = 1/3. 10 

10 In order to correctly reproduce the dynamics for small three-volume, namely, the stalk 
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Reverting back to a formulation in terms of the scale factor a(r), and by 
suitable rescaling of r and o(t)), we can now write the effective action of causal 
dynamical triangulations at fixed continuum four-volume V4 as 



<eff _ 




) 



(40) 



Again, this expression is valid for sufficiently large scale factor, since we have 
explicitly dropped any reference to the quantum corrections at small a. Compar- 
ing now with (jSHJ), the startling conclusion is that we have indeed rederived the 
Euclidean minisuperspace action. However, here we have done it from first prin- 
ciples, and up to an overall sign. The strongest evidence so far that the effective 
action is indeed given by ()40|) comes from using it as input to create an artificial 
distribution of three- volumes N 3 (t) [2U|. We then evaluate the volume- volume 
correlator c ^ (x) of eq. (jl(Jj) on this ensemble in the same way as we would do for 
genuine Monte Carlo data. As illustrated in Fig. the matching of this new 
curve with the previous Monte Carlo results (c.f. Fig. |7J) is impressive. 

One immediate conclusion is that the collective effect of the dynamical de- 
grees of freedom other than a(r), which are ignored in minisuperspace models, is 
to change the sign of the kinetic term for the scale factor (the "global conformal 
mode") from —aa 2 to +aa 2 . This is the same nonperturbative cure of the confor- 
mal divergence of the Euclidean path integral already encountered in dimension 
3 [33J EH] and argued for from a continuum point of view in [56] . It provides 
further evidence that the well-known sicknesses and ambiguities of the Euclidean 
path integral approach are features of the semiclassical and symmetry-reduced 
approximations, not shared by a full nonperturbative treatment which is rooted 
in a causal, Lorentzian formulation. 

We have begun an exploration of the consequences of our result for cosmology 
and quantum cosmology. When searching for a physical interpretation, one must 
keep in mind that the effective action (|4(J|) we extracted directly from the com- 
puter simulations holds for fixed four-volume. Translating between propagators 
of fixed four-volume and fixed cosmological constant requires a further Laplace 
transformation [20] • Moreover, we are still in the Euclidean sector of the theory, 
and must Wick-rotate back to obtain physical, Lorentzian results. 

Nevertheless, it is exciting to see that the computation of some physical quan- 
tities is already within reach: for example, we showed in [20] how to derive the 
semiclassical wave function of the universe without explicitly performing an in- 
verse Wick rotation. The wave function can be compared directly with those 

observed at large times r, the function N^ 3 in l|38|) with a = 1/3 must be replaced by a 
function of N3 whose derivative at goes like N% , v > 0, as explained in [20|. A simple choice 
is the replacement iV 3 1/3 ->• (1 + iV 3 ) 1/3 - 1. 
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Figure 13: The scaling of the volume- volume distribution, as function of the rescaled 
time variable x = ^/(-/V^) 1 / 4 , taken from Fig. [7[ but now superimposed with the cor- 
relator obtained when using the minisuperspace action (|40|) to generate a distribution 
of three-geometries (black curve) . 

arising in standard minisuperspace treatments, like those of Vilenkin |4T[ 
Hawking ^SJ HH] and Linde [3B] • To relate to the dynamics of the real universe, 
especially at very early times, we will need to couple matter to our causal geome- 
tries. These issues are currently being investigated and will be reported in the 
near future. 



. 8 
. 6 

^? 

u0 . 4 
0.2 





6 Nonclassical structure of the universe: 
The geometry of spatial slices 

Having established that our quantum universe possesses certain classical features 
at sufficiently large scales, we are interested in a more detailed characterization 
of its geometry. We report in this section on various measurements associated 
with the geometry of spatial slices r = const, and in the next section on that of 
"thick spatial slices", by which we mean geometries of time extension At = 1, 
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which are "sandwiched" between adjacent spatial slices of integer time r. Our 
measurements describe certain invariant geometric properties of the ground state 
of spacetime, but we make at this stage no attempt to relate them to actual 
observables that one could go out and measure in the real universe. They will 
merely serve to illustrate the complexity of the geometric structure and highlight 
some of its nonclassical features. 

Because our universes have a well-defined global notion of (proper) time r, it 
is relatively straightforward to perform measurements within a slice of constant 
r. We first restrict ourselves to slices S of integer-r. Such slices consist entirely 
of spatial tetrahedra (whose edges are all spacelike). The slices are of the form 
of spatial interfaces of topology S* 3 between adjacent sandwiches made out of 
four-simplices. As we will see, the Hausdorff dimension of the slices turns out 
to be three, as one would have naively expected in a four-dimensional universe. 
However, the slice does not behave three-dimensionally under diffusion, and we 
will see that its structure is fractal in a precise sense which will be described 
below. 



6.1 Hausdorff dimension 

Our algorithm for determining the Hausdorff dimension d h of the spatial slices is 
as follows. From a given S of discrete volume N 3 , we randomly pick a tetrahedron 
i . From i , we move out by one step and obtain n(l) tetrahedra at distance 1 
from i (that is, its nearest neighbours). Moving out by a further step, there will 
be n(2) tetrahedra at distance 2, and so on until all tetrahedra have been visited. 
The numbers n(r) recorded for each distance r sum up to 

fmax 

E < r ) = N *- ( 41 ) 

Finally, we measure the average linear extension 

3 r 

of the spatial slice. For a slice of volume N 3 , this process is repeated iV" 3 /50 + 1 
times with different randomly chosen initial tetrahedra. The expectation value 
(r) is then averaged over these measurements. In this way, we obtain for every 
slice a data pair 

{{{r)),N 3 }, (43) 

representing one "bare" measurement. 

This process is performed for all t spatial slices of the geometry, thus collecting 
t data pairs. We have typically done this for about 1000 different spacetime 



30 



geometries, obtaining altogether lOOOt measurement points. The final results are 
sorted by their ((r))-value (in practice a continuous variable) and averaged over 
a sequence of 100 consecutive points. This reduces the number of data points 
to lOt, which are then displayed on a log-log plot. In the presence of finite-size 
scaling, we expect them to follow a curve 

(N 3 )(r)<x((r)+r ) d \ (44) 

defining a spatial Hausdorff dimension dh- In relation (jUj) , we have included 
a finite shift r of the linear size r, based on earlier experience with finite-size 
corrections in the scaling relations |3TH 14*3*] . Our results, which show universal 
behaviour and scaling, are presented in Fig. HU The shift value r = 1.75 gives 
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Figure 14: Log-log plot of the average linear geodesic size (r) jy 3 versus the three- volume 
N%, measured for A = 0.4. The straight line corresponds to a Hausdorff dimension 
dh = 3. Similar measurements for A = 0.5 and 0.6 yield virtually indistiguishable 
results. 

the best linear fit. The resulting value of dh ~ 3 is practically independent of tq. 
6.2 Spectral dimension 

Given the results of the previous sections, one might be tempted to conclude 
that the geometry of our dynamically generated ground state simply is that of 
a smooth four-dimensional classical spacetime, up to Gaussian fluctuations. A 
more detailed analysis of the geometry of spatial slices makes explicit that this 
is not so. - We start by measuring the spectral dimension d s of the spatial slices 
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by studying diffusion on them. Analogous to the four- dimensional case described 
in Sec. 15.21 we choose an initial probability distribution 



K T (i,i ;a = 0) = 5,, 



(45) 



on the three-dimensional triangulation T of the slice S, which is zero everywhere 
except at a randomly chosen tetrahedron io, and define the evolution in discrete 
time a by the rule 



where k — > j labels the (four) nearest neighbours of the tetrahedron j. We 
measure the return probability Prii^cr) = ifr(io, zo; &)■ The measurements of 
the return probabilities are repeated a number of times for each slice S with 
different starting points zq, precisely as in the measurement of the Hausdorff 
dimension of S. They are then repeated for independent configurations (four- 
dimensional triangulations generated by the Monte Carlo- procedure), again with 
the same statistics as in the measurement of the Hausdorff dimension of the three- 
dimensional spatial slices. Results are stored for averaged slices with volumes 
N 3 = n x 500 ± 50, because we expect them to depend on the system volume. In 
this way we define the average return probability P Ns (a) and extract the spectral 
dimension 



from the logarithmic derivative exactly as in the four- dimensional diffusion pro- 
cess. 

The result is displayed in Fig. for even and odd times separately, because 
of the different behaviour of P^ 3 (a) for short evolution times. As can be seen 
the result is very stable for larger diffusion times and gives d s = 1.56 ± 0.1 for the 
spectral dimension, which is only about half the expected "classical" value! 

6.3 Critical exponents 

The noncanonical value of the spatial spectral dimension shows that the quantum 
geometry has a nontrivial microstructure, which definitely is not that of a smooth 
four- dimensional manifold. In fact, we already have come across further evidence 
of this fact in the measurement of the spacetime spectral dimension which ap- 
proached the value four only at long distances. One way to obtain an additional 
characterization of geometries at a fixed time is through the functional form of 
their "entropy" or "number of states" j\f(N) as function of their discrete volume 




(46) 




(47) 
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Figure 15: The logarithmic derivative of the return probability for kq = 2.2 and A = 0.4, 
and spatial volumes N% between 500 and 2000. The measurement was performed for a 
system with N4 = 91.1k and £ = 40. The increasing and decreasing lines correspond to 
even and odd times, and illustrate the different behaviour for short diffusion times a. 



N . Two different functional forms for the entropy have appeared in previous 
studies of simplicial three- dimensional quantum gravity [4*4*] . namely, 

M(N 3 ) oc jV 3 - 3+ ^e wAr3 , (48) 
Af(N 3 ) oc e ^ o7V3 - CliV 3 + -, v < 1, (49) 

with a nonuniversal constant hq measuring the leading exponential growth. The 
critical exponent 7 str in f|4*8*j) characterizes the fractal structure of the geometry 

0"""] |("U] . The larger the value of 7 str , the more pronounced the tendency for a 
given three-geometry to develop so-called baby universes 11 . 

The functional form (|4*9*|) on the other hand indicates a suppression of baby 
universes and formally may be seen as the limit 7 str — ► — oo of (J48|) . Below we will 
meet both functional forms. It is important to understand that the "entropy" 
Af(N 3 ) in (*4*8*j) and ("4*9*|) does not merely count the possible three-geometries, but 
also includes the weight coming from the action: in counting geometries we are 

11 In the present context, a "baby universe" is a part of the three-dimensional spatial tri- 
angulation which is connected to the remainder by a "minimal neck" . This means that the 
interface between the triangulated baby universe and the "parent universe" is topologically a 
sphere and of minimal size. Equivalently, it is given by the surface of a tetrahedron consisting 
of four triangles. This is easier to visualize in a two-dimensional simplicial manifold where the 
analogous minimal neck is an S 1 consisting of three edges. 
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making an importance sampling where the action defines the probability weight 
of the geometry. We saw above that the choice of bare coupling constants in 
the four-dimensional theory determines in which of the different phases of the 
lattice theory we end up. The "entropy" of the spatial geometries at fixed time 
is different in the three phases, even if the underlying geometries are identical, 
and simply given by the space of all (abstract) triangulations of S 3 . 

If the functional form of the entropy is given by eq. ()48j1 . a convenient way of 
measuring the exponent 7 str is through the distribution of baby universes (23 EDI • 
Calling B the volume of the baby universe (the number of tetrahedra contained 
in it) and N 3 the total volume of the spatial geometry, we expect a distribution 
of the form $R ED] 

P(x) oc x~ 2+ ^ r (l + + •••], (50) 



where 



X = ^( 1_ ^)' ° <X<1 ^ W 

and c(N 3 ) oc 1/N 3 is a volume-dependent constant. Relation (JHOJl is expected 
to hold for sufficiently large B. (It was derived under the assumption that the 
samples of both baby and parent universes have an asymptotic distribution of the 
form (jUJ.) Since it is very difficult in the present set-up to collect data for any 
particular fixed value of N 3 , we have performed "integrated" measurements to 
determine 7 str , by collecting the distribution of the number of baby universes of a 
given x for various volumes A3 > 400. Also in this case we expect the distribution 
of baby universes to follow 

P(x) oc x~ 2+lstt ^1 + ^ + . . ?j , (52) 

with a constant d depending now on the range of three- volumes A3 over which we 
have averaged the distribution. The resulting distributions for k = 2.2, A = 0.4 
and total system sizes A4 = 45.5, 91 and 181k are shown in Fig. Although 
the measurement is of rather poor quality, and - depending on which range of x 
is used for the fit - the uncertainty in 7 str rather big, the approximate linearity of 
the plots strongly suggests that the distribution is indeed dominated by a power 
law for x. Fitting 7 str to formula (p^j) . we obtain 7 str = 0.35 ± .09. 

Another, seemingly independent, way of characterizing the short-distance be- 
haviour of the spatial slices S is by measuring the coordination number of vertices 
in S, i.e. the number n v of spatial tetrahedra sharing a given vertex v. 12 In the 



12 Note that for a 3d manifold there is a linear relation between the coordination number of 
a vertex v and the number of links emerging from v. 
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Figure 16: The distribution of the number of baby universes with minimal neck as 
a function of their rescaled size x, for kq = 2.2 and A = 0.4. Distributions are nor- 
malized by the total number of baby universes in constant-time slices with N3 > 400. 
This number grows with the total spacetime volume N4, which is 45.5k (bottom), 91k 
(middle) and 181k (top curve). 



simulations, we observe vertices with relatively high coordination numbers n, and 
a distribution of n which can be approximated a simple power law 

P{n) oc —— . (53) 

As usual, we have included a finite shift uq to take into account finite-size effects. 
For a three-dimensional triangulation, n > 4 and n may take only even values. 

Fig. El illustrates the distribution of coordination numbers for 4 < n < 200, 
measured for a system with Ko = 2.2, A = 0.4, t = 40 and system size N4 =181k. 
The data points have been collected from spatial slices with three-volume N3 > 
500, and can be seen to behave perfectly linearly on the log-log plot. The straight 
line in Fig. El corresponds to a fit with a = 1.53. In theories of random geometry, 
this type of power-law behaviour of the vertex coordination number was first 
observed in dynamically triangulated bosonic string theory [S3] . The value of the 
exponent a extracted from the data by x 2 -flts depends slightly on the choice of 
no and of the cut-off in the data. Our estimate is a = 1.53 ± 0.05. 

Another way to determine a is by studying the distribution of points with the 
largest coordination number as a function of the spatial volume A^ 3 (or the total 
number of vertices NqIN^), the two numbers being proportional for the phase of 
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Figure 17: Distribution P(n) of the coordination numbers n of vertices in a three- 
dimensional spatial slice, for kq = 2.2 and A = 0.4. The straight line represents a 
numerical fit to formula with no = 2. 

four- dimensional gravity we are considering here). One expects a behaviour of 
the form 

P{n > n max (T)) oc (54) 

where n max (T) denotes the maximal coordination number for a particular three- 
dimensional triangulation T. The proportionality ()54j) is expected to be valid 
if the distribution P{n) is such that a typical n max {T) is so small that it is 
not influenced by the fact that the volume N 3 (T) contains only a finite number 
N (T) of vertices. (Clearly, n max (T) cannot be larger than N (T)/2.) Under 
this assumption the orders n v of the N (T) vertices v in T are representative of 
the distribution P(n). Since there are N Q (T) of them, the chance of "drawing" 
an (abstract) vertex from the distribution P(n) with a coordination number n > 
n max (T) is of the order 1/N (T), which, given ^(-^3) oc N 3 , is equivalent to (JHU). 
Consequently, using the scaling (JoHj) . we obtain that 

rw(iV 3 ) oc Nl ,a . (55) 

To check this property, we measure for each surface S the largest coordination 
numbers in S. Since these fluctuate considerably, we measure the average of the 
three largest values for each S. We then follow a strategy similar to that adopted 
for the measurement of the Hausdorff dimension. Namely, we collect data pairs 
{N 3 ,n max }, sort the resulting set by volume N 3 and average over a set of 1000 
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measurements. Fig. IT81 shows the result in the form of a log-log plot, confirming 
the expected linear behaviour. The solid line corresponds to the value 
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Figure 18: Distribution of the largest coordination number n max as a function of the 
spatial volume N 3 , for k = 2.2, A = 0.4, t = 40 and iV 4 = 91.1A;. 



in fair agreement with the previous measurement from smaller correlation num- 
bers, and the fit gives 1/a = 0.681 ± .05, translating roughly into a = 1.468 ± 0.12. 
The origin of the uncertainty is the chosen range in the three-volume N 3 taken 
into account in the fitting. 

6.4 Relation to branched polymers 

The above measurements lead to the intriguing conclusion that all the geometric 
properties of the spatial slices measured so far can be modelled by a particular 
kind of branched polymers 13 , as we will now go on to explain. 

A (planar) branched polymer is an (abstract) one- dimensional polymer which 
admits arbitrary branching. It is modelled by (abstract) connected planar graphs 
which contain no loops. As a consequence, cutting any link of such a graph will 
result in two disconnected graphs. At each vertex an "incoming" link is allowed 

13 In combinatorics, branched polymers are conventionally called rooted trees, where "rooted" 
means that one link or vertex is marked, and therefore distinguished from the rest. 



1 



(56) 



0.681 



a 
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to "branch" into n— 1 "outgoing" links with probability p n . The statistical model 
of (rooted) branched polymers is then defined by the partition function 

Zbp(v) = ^e-^^flp(JV) (57) 

N 

N 

Zbp(N) = ^p(BP N ), p(BP N ) = l[p n , (58) 



BP N re=l 



where the summation is over all branched polymers with N vertices, and the 
weight p{BP) is the product of the branching probabilities. One can also embed 
the abstract branched polymers in a target space, either a lattice or R d . In the 
latter case, it is convenient to work with a quadratic action (x u — x v ) 2 for pairs of 
nearest-neighbour vertices u and v with coordinates x u and x v , and perform an 
additional integration over the target-space coordinates in the partition function. 

The fractal properties of branched polymers are well understood. The leading 
behaviour for large A" is an analogue 14 of eq. (|48|) . namely, 



Z BP {N) ~ iV 7 " 2 e»° N , (59) 
where 7 can take a range of real values, depending on the branching probability 



p n . The intrinsic Hausdorff dimension <il of branched polymers with 7 > is 



E21 E3] 

a h — 

7 



4 n) = ~, (60) 



with spectral dimension given by [501 

ds = 7^-- (61) 
1 + 7 

We note that formulas (j6Up and ()61|) already provide a close match with our 
data if we identify 7 with the 7 str of the simulations. As we will see, there 
is indeed a concrete branched-polymer model which has the 7-value we have 
observed. In regularized bosonic string theory one is naturally led to studying 
branched polymers with branching weight [64J 

p n ~ — for n > I. (62) 

Generically, the value of 7 for a given a is 

7(a) = 1/2 for a < 1, , g3 * 
7(a) = 1 — a for 1 < a, 



14 The difference between the power 7—2 in eq. I|59|l and 7—3 of eq. (|48|l arises because of the 
use of rooted branched polymers, which gives rise to an extra factor of N. 
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which does not cover the 7-range we are interested in. However, if the partition 
function satisfies a certain technical condition 15 , the relations (|63|) are changed 
to 

7(a) = (a— l)/a for l<a<2 . . 

7(a) = 1/2 for 2<a ^ } 

and see [TH] for a detailed discussion). While the associated condition (IM1) 
seems rather artificial from the perspective of branched polymers, it appears in 
a natural way if the branched polymers emerge as a limit of higher-dimensional 
models like noncritical strings with central charge c > 1. As a matter of fact, the 
generic value of 7 in this case is 1/3 [H2]- Summarizing the critical exponents for 
this particular class of branched polymers, 



i=\i 4 m) = 3 > ds = % a = \> ((iC) 



we see that within the measuring accuracy, they are precisely the numbers char- 
acterizing the spatial geometries of causal dynamical triangulations at constant 
time. Although the vertex order n v of our three-dimensional triangulations has 
a priori nothing to do with the branching of polymers, the above results are a 
strong hint that there effectively is such a connection. While this confirms our 
earlier assertion that the geometry of the spatial slices behaves very nonclassi- 
cally, it opens the exciting possibility that quantum gravity as described by this 
model may be amenable to analytic approximations, precisely as was the case in 
the study of noncritical strings with c > 1 



7 Nonclassical structure of the universe: 
The geometry of "thick" spatial slices 

Another substructure of spacetime whose geometry can be investigated straight- 
forwardly are what we shall call "thick" spatial slices, namely, the sandwiches 
of geometry contained in between two adjacent spatial slices of integer times r 
and r + 1, with topology I x S 3 . As illustrated in Fig. [TJ they are made up of 
four-simplices of types (4,1) and (3,2) and their time-reversed counterparts (1,4) 
and (2,3). Because of their finite time extension, one might expect the geometry 
of the thick slices to behave more "classically", and this is indeed corroborated 



15 For (|(j5fl to hold one needs 

= 00, (64) 



dZ BP {fi) dZ B p{n) 
< tor /! > /i , 



dfj, dfi 



which is naturally satisfied for some branched polymer models coming from noncritical string 
theory [H^]. 
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by our simulations, which include a measurement of the Hausdorff and spectral 
dimension of the thick slices. 

A convenient way to represent a sandwich geometry is by considering the 
three-dimensional piecewise flat geometry that results by cutting the thick slice 
at time r + 1/2. The analogous procedure in one dimension less, which is eas- 
ier to visualize, gives a generalized two-dimensional triangulation consisting of 
triangles and squares [SH ESI E2] The triangles are the intersection surfaces 
of (3,1)- and (l,3)-tetrahedra, and the squares those of (2,2)-tetrahedra. In 
one dimension higher, the fundamental building blocks are no longer tetrahedra 
but four-simplices. Consequently, the intersection patterns are three-dimensional 
flat building blocks, namely, spatial tetrahedra from slicing the (4,1)- and (1,4)- 
simplices, and spatial triangular prisms from slicing the (3,2)- and (2,3)-simplices. 
In order to distinguish between a tetrahedron that comes from cutting a (4,1)- 
simplex from one associated with an upside-down, (l,4)-simplex, one colour-codes 
the links of the intersection hypersurface into red and blue, say. The red links are 
those that will be "shrunk away" when the constant-time surface approaches the 
integer-time r+1, and similarly the blue links will disappear as the hypersurface 
is moved towards time r. As a result, the tetrahedra have all links of one colour, 
whereas the prisms are always bi-coloured (see also [221 )■ Note that the average 
spatial volume in terms of counting three-dimensional building blocks is larger 
(by a factor of about 2.5) at slices of half-integer time, because - unlike at integer 
times - all four-simplices between r and r + 1 contribute to the counting. 

7.1 The Hausdorff dimension 

To determine the Hausdorff dimension du of the thick slices, we pick a ran- 
dom four-simplex from such a slice, or, equivalently, a random three-dimensional 
building block from the corresponding intersection pattern at half-integer time. 
Continuing in this three-dimensional representation of the slice geometry, we pro- 
ceed exactly as we did when measuring the Hausdorff dimension d^ of the spatial 
slices at integer-r. Starting from the randomly chosen initial three-dimensional 
building block, we move out, one step at a time, to the nearest neighbours at 
distance r, from the previous set of building blocks at distance r — 1. In terms of 
four-simplices, this means we never move out to simplices earlier than r or later 
than r+1, that is, by definition we never leave the thick slice. At each step, we 
keep track of the number n(r) of building blocks, whose sum is equal to the total 
slice volume N, 




(67) 



40 



and which enables us to compute the average linear size 

( r ) = Jj^2 rn ( r ) (68) 

r 

of the thick slice. For a given slice of volume N, this process is then repeated 
N4/IOO + 1 times with randomly chosen initial building blocks. The average size 
(r) is again averaged over this set and the data pair {({r)),N} stored. Using 
the same procedure as described following formula (J4*H|) above, we extract the 
Hausdorff dimension du from the relation 

N A oc ((r) +r ) dH } (69) 

with a finite constant r taking into account finite-size corrections. The results are 
presented in Fig. E3in the form of a log-log plot. The straight line corresponds 
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Figure 19: Log-log plot of the average linear geodesic size (r)jv versus the volume N 
of a thick slice. The straight line corresponds to a Hausdorff dimension d# = 4. 

to du = 4, and the best fitted value is du = 4.01 ± 0.05. Measurements were 
performed on spacetimes with total size N4 =20, 40 and 80k, at Kq — 2.2, A = 0.6 
and t = 80. - We conclude that the transition from a genuine constant-time 
slice to a thick slice adds one extra dimension to the Hausdorff dimension - the 
thickened slice already "feels" some aspects of the four-dimensionality of the full 
spacetime. Note that this would not be true for a thick slice of a classical, regular 
lattice, whose Hausdorff dimension would be the same as that of a "thin" slice. 
The likely mechanism for the dimensional increase in the triangulations is the 
appearance of "short-cuts" between two tetrahedral building blocks, say, once 
one moves slightly away from a slice of integer-r. 
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7.2 The spectral dimension 

Next, we turn to the spectral dimension ds of the thick slices. When adapting 
our earlier measurement of the spectral dimension d s of spatial slices, we have to 
take into account that the tetrahedral building blocks in the three-dimensional 
intersection pattern at half-integer r have four neighbours, whereas the prisms 
have five. We start from an initial probability distribution 

K T (i,i ,a = 0) = 5 i>io (70) 

concentrated at a randomly chosen three-dimensional building block with label 
io in the three-dimensional "triangulation" T (consisting of tetrahedra and trian- 
gular prisms), and define the evolution process by 

K T (j,i ;a + l) = V— K T (k,i Q ;a), (71) 

k^-j 

where denotes the number of neighbours of building block i, and k — > j are 
the nearest neighbours of the "simplex" j (which can now also be a prism). By 
construction, the total probability 

Y^K T (i,i ;a) = l (72) 

i 

is conserved. As before, we measure the return probability KT(io,io\cr). After 
the appropriate averaging over initial points z'o and triangulations T we obtain 
the quantum return probability Pn{o) as a function of the number of building 
blocks N in the triangulations T (we keep the number iV approximately constant 
in the triangulations included in the average). The expected behaviour 

Pn(ct) oc a~ ds ^ 2 + finite-size corrections (73) 

defines a spectral dimension ds, which in general will be different from the spectral 
dimension d s determined for the constant-time slices. The value of ds is extracted 
by measuring the logarithmic derivative of Pjv(c) as a function of the evolution 
time a, 

is(<j) = -2 f^M. ( 74 ) 
a log a 

The results of our measurements are displayed separately for even and odd times 
in Fig. |20J since the behaviour of Pat (a) is different for very small a (as already 
observed in our previous measurements of spectral dimensions). As seen in the 
figure, beyond small a we do not find the perfect constant logarithmic derivative 
observed earlier for the "thin" slices at integer time. If we nevertheless treat 
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Figure 20: The logarithmic derivative of the return probability for ko = 2.2 and A = 0.6 
and thick slices of volumes N between 900 and 2000. The measurement was performed 
for a system with iV~4 = i61&; and i = 80. Curves extrapolating between even evolution 
times decrease rapidly at small a, whereas curves of odd times stay above ds = 2 at 
small times. 

the output curves as constants and take their mean value as representative of the 
logarithmic derivative ds (assuming it is cr-independent), we obtain ds = 2. 0±0. 2. 
Although it lies above the value d s for the spectral dimension of the spatial slices, 
this dimensionality is still highly nonclassical, and reflects the presence of some 
intermediate regime between a purely spatial and a purely spacetime geometry. 

7.3 Baby universes 

To further characterize the geometry of thick slices, we also looked into its baby 
universe structure. Minimal necks can take the form of either four triangles, form- 
ing the surface of a tetrahedron, or of two triangles and three squares, forming the 
surface of a triangular prism (of course, without an interior tetrahedron or prism 
being located there). It turns out that the four- dimensional manifold constraints, 
which ensure that the four-simplices are glued in a way that results in a simplicial 
manifold, impose severe restrictions on the types of baby universes with minimal 
necks that can occur. 16 For a tetrahedral minimal neck, the baby universe which 
branches off there can only contain building blocks which are tetrahedra of the 

16 Analogous manifold conditions for intersection graphs in three-dimensional simplicial space- 
times were formulated in 33 ■ 



43 



same colour. Minimal necks of the form of a prism cannot occur altogether. We 
attempted to measure the distribution of "monochrome" tetrahedral baby uni- 
verses, but found them to be small in size, within a range practically independent 
of the slice volume N. We did not check for baby universes with larger necks. 
Our tentative conclusion is that the microscopic structure of the thick slices is 
more regular than that of the genuine spatial slices of fixed constant time. 

It is tempting to conjecture that the absence of baby universes means that 
the entropy function M{V) for the thick spatial slices has the functional form 
()49|) rather than (|48p. which was the form observed for the thin spatial slices at 
fixed integer time. As mentioned in Sec. I6.3| the entropy behaviour (J4*§|) signals 
an almost complete suppression of baby universes. Moreoever, if the exponent 
v in eq. (f4*HJ) was equal to 1/3, it would fit neatly with the observation that the 
effective macroscopic action (|4()J). which reproduces so well the macroscopic four- 
dimensional features we observe in the simulations, does contain the characteristic 
linear term in the scale factor a(r) (corresponding to a factor V^ 3 (t) in terms of 
the continuum three- volume). It would also be in line with our expectation that 
the structure of the thick spatial slices, rather than that of the thin slices at fixed 
time should relate most directly to the global four-dimensional structure. 

8 Summary and outlook 

This paper describes the currently known geometric properties of the quantum 
universe generated by the method of causal dynamical triangulations, as well as 
the general phase structure of the underlying statistical model of four-dimensional 
random geometries. The main results are as follows. An extended quantum 
universe exists in one of the three observed phases of the model, which occurs for 
sufficiently large values of the bare Newton's constant G and of the asymmetry 
A, which quantifies the finite relative length scale between the time and spatial 
directions. In the two other observed phases, the universe disintegrates into a 
rapid succession of spatial slices of vanishing and nonvanishing spatial volume 
(small G), or collapses in the time direction to a universe that only exists for an 
infinitesimal moment in time (large G, vanishing or small A). In either of these 
two cases, no macroscopically extended spacetime geometry is obtained. 

By measuring the (Euclidean) geometry of the dynamically generated quan- 
tum spacetime 17 in the remaining phase, in which the universe appears to be 
extended in space and time, we collected strong evidence that it behaves as a 
/ewr-dimensional quantity on large scales. First, finite-size scaling of the volume- 
volume correlator for spatial slices is observed when the time variable is rescaled 

17 Since our universe is a weighted quantum superposition of geometries, all "measurements" 
refer to expectation values of geometric operators in the quantum ground state. 
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by the fourth root of the discrete four- volume, eq.(|ll|). within measuring accu- 
racy. Second, the large-scale spectral dimension of the universe, determined from 
a diffusion process on the quantum ensemble of geometries, is compatible with 
the value 4. 

The dynamical derivation of the correct classical result for the dimension from 
a completely nonperturbative path-integral setting provides a considerable boost 
to our efforts to construct a theory of quantum gravity by the method of causal 
dynamical triangulations. Furthermore, we were able to show that the universe 
possesses other features of classical geometry at sufficiently large scales. Namely, 
using the four-dimensionality at large scales, we derived an effective Euclidean 
action for the scale factor of the universe (equivalently, the volume of its spatial 
slices). For the case of fixed spacetime volume, this action coincides up to an 
overall sign with a simple minisuperspace action frequently used in quantum 
cosmology. Again, we believe that this is a highly nontrivial and very exciting 
result. It has already enabled us to derive a "wave function of the universe" from 
first principles [20J. The change in sign of the action implies that the Euclidean 
sector of our quantum gravity model is well-defined. This is a welcome result 
since, unlike in standard quantum cosmology, no ad hoc cure is needed to deal 
with a conformal kinetic term that renders the action unbounded below. Some 
care is required in interpreting it, since we still need to perform an inverse Wick 
rotation to get back to the physical, Lorentzian sector of the theory. A detailed 
discussion will appear elsewhere. 

Returning to the analysis of the quantum geometry of our universe, we gained 
further insights by performing measurements of a more local nature by looking 
at the intrinsic geometry of "thin" and "thick" spatial slices. A thin slice is by 
definition a spatial geometry at fixed integer time somewhere inside a universe, 
and a thick slice a piece of spacetime contained between two adjacent thin slices. 
We found that in general the measured observables exhibit nonclassical behaviour, 
in the sense that dynamically determined dimensionalities of the slices do not 
coincide with the values one would have expected naively for the corresponding 
submanifolds of a classical four-dimensional spacetime. This result is not in 
contradiction with the four-dimensional behaviour of the universe at large scales, 
but underscores the fact that the local geometry of the quantum spacetime has 
highly nonclassical features. 18 

We found that the measured Hausdorff dimension of thin slices comes out 
as 3 with good accuracy, but that the spectral dimension is only around 1.5! 
By making additional measurements of the distributions of vertex coordination 
numbers and of baby universes in the thin slices we were able to show that 

18 One should keep in mind that slices of constant time are in no way physically distinguished 
in a theory of pure gravity and their geometric properties therefore not directly related to 
physical observables. 
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the entire set of parameters is matched by the critical exponents of a particular 
class of branched polymers, potentially opening a new window to an analytic 
investigation of quantum spacetime. 

The thick slices were found to have Hausdorff dimension 4 within measuring 
accuracy, and a spectral dimension of around 2. This, and the apparent absence 
of baby universes of appreciable size indicate that the geometry of the thick slices 
represents a genuinely intermediate case between that of the thin slices and of the 
full spacetime. In particular, their geometry - although far from being classical 
- seems to be more regular than that of the thin slices. 

The most local measurement of quantum geometry so far is that of the spectral 
dimension of spacetime at short distances, which provides another quantitative 
measure of the nonclassicality of geometry. As we have seen, the spectral dimen- 
sion changes smoothly from about 4 on large scales to about 2 on small scales. Not 
only does this (to our knowledge) constitute the first dynamical derivation of a 
scale-dependent dimension in full quantum gravity, but it may also provide a nat- 
ural short- distance cut-off by which the nonperturbative formulation evades the 
ultra-violet infinities of perturbative quantum gravity. Remarkably, evidence for 
a dimensional reduction from four to two has been found recently in the so-called 
exact renormalization group approach to (Euclidean) quantum gravity, truncated 
to a finite-dimensional parameter space 2]. One could be particularly encouraged 
by the fact that despite very different starting points and quantization methods 
the two approaches nevertheless seem to lead to compatible results. 

In summary, what emerges from our formulation of nonperturbative quantum 
gravity as a continuum limit of causal dynamical triangulations is a compelling 
and rather concrete geometric picture of quantum spacetime. Quantum space- 
time possesses a number of large-scale properties expected of a four-dimensional 
classical universe, but at the same time exhibits a nonclassical and nonsmooth 
behaviour microscopically, due to large quantum fluctuations of the geometry at 
small scales. These fluctuations "conspire" to create a quantum geometry that is 
effectively two-dimensional at short distances. 

Of course, our large-scale measurements of geometry so far do not characterize 
the classical properties of our quantum spacetime completely, and further work 
is under way to determine its large-scale curvature properties, and check that 
Newton's inverse square law can be recovered in an appropriate limit (see 
for related previous attempts). Having established the existence of a meaningful 
classical limit opens the door to the ultimate subject of our interest, the quan- 
tum behaviour of the local gravitational degrees of freedom. On the one hand, 
together with appropriate matter degrees of freedom, they should be able to de- 
scribe the quantum behaviour of the very early universe, which may have left an 
imprint on the cosmic microwave background radiation we observe today. On the 
other hand, they will characterize the fluctuations around a large-scale classical 



46 



background, like that of our present universe, leading to quantifiable deviations 
from the predictions of classical relativity. 
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